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Context. Recent work on several /? Cephei stars has succeeded in constraining both their interior rotation profile and their convective 
core overshoot. In particular, a recent study focusing on Ophiuchi has shown that a convective core overshoot parameter of a m , = 
0.44 is required to model the observed pulsation frequencies, significantly higher than for other stars of this type. 
Aims. We investigate the effects of rotation and overshoot in early type main sequence pulsators, such as ft Cephei stars, and attempt 
to use the low order pulsation frequencies to constrain these parameters. This will be applied to a few test models and the /? Cephei 
star 9 Ophiuchi. 

Methods. We use the 2D stellar evolution code ROTORC and the 2D linear adiabatic pulsation code NRO to calculate pulsation 
frequencies for 9.5 M Q models evolved to an age of 15.6 Myr. We calculate low order p-modes (I < 2) for models with a range of 
rotation rates and convective core overshoot parameters. These low order modes are the same range of modes observed in 9 Ophiuchi. 
Results. Using these models, we find that the convective core overshoot has a larger effect on the pulsation frequencies than the rota- 
tion, except in the most rapidly rotating models considered. When the differences in radii are accounted for by scaling the frequencies 
by V(GM/i?(40°) 3 ), the effects of rotation diminish, but are not entirely accounted for. Thus, this scaling emphasizes the differences 
produced by changing the convective core overshoot. We find that increasing the convective core overshoot decreases the large sepa- 
ration, while producing a slight increase in the small separations. We created a model frequency grid which spanned several rotation 
rates and convective core overshoot values. We used this grid to define a modified,^ 2 statistic in order to determine the best fitting 
parameters from a set of observed frequencies. Using this statistic, we are able to recover the rotation velocity and convective core 
overshoot for a few test models. We have also performed a "hare and hound" exercise to see how well ID models can recover these 
parameters. Finally, we discuss the case of the fj Cephei star 9 Oph. Using the observed frequencies and a fixed mass and metallicity, 
we find a lower overshoot than previously determined, with a m , = 0.28 ± 0.05. Our determination of the rotation rate agrees well with 
both previous work and observations, around 30 km s -1 . 

Key words, stars: early-type - Stars: oscillations - Stars:individual (9 Oph) - Stars: rotation - Stars: interiors 



1. Introduction 

Recently, great progress in the asteroseismology of /3 Cephei 
stars has been made thanks to extensive observational cam- 
paigns, which have allowed constraints to be place d on the in- 
terior p rop erties of some of the se stars. For example, Aerts et al. 
d2003l) and I Aerts etaf] d2004bh have compiled and analyzed 21 
years of photometry for V836 Cen, identifying six frequencies 
and their degree and order. Subsequent modeling has placed con- 
straints on the mass, age, metallicity, convective co re overshoot- 
ing (a m ) and internal rotation profile for this star (Dupret et al. 
2004). They find strong evidence for a ov — 0.1 + 0.05 in the 
absence of rotational mixing. Although the rotation rate for this 
star is quite slow, around 2 km s , the observed frequencies can 
not be matched with a uniformly rotating model. 

A second Cephei star, v Eridani, has also been studied ex- 
tensively with both photometric and s pectroscopic campaigns 
dHandler et all 12004: Ae rts et al.ll2004ah . Nine modes were de- 
tected for this star, including th e radial mode and two £ = 
1 triplets (De Ridd er et alJ I2004T) . Modeling of v Eridani has 
also shown that non-un i formly rotating mode ls are required 
dPamvatnvkh et al.l 12004 lAusseloos et alj|2004l) . As for V836 
Cen, this star cannot be uniformly rotating to match the observed 
frequencies, and some convective core overshooting may be re- 
quired. There is also some indication that the interior chemical 



composition is not homogeneous, with Fe overabundant in the 
driving zone. 

Several other [5 Cephei stars have also been s uccessfully 
model ed. Using observations by the MOST satellite. [Aerts et all 
(2006) were able to place constraints on the physical parameters 
of 6 Ceti, including constraining the convective core overshoot 
(a nv = 0.2 + 0.05). Seismic modeling of /3 CMa has constrained 
the core overshoot to a ov = 0.2 ± 0.05, as well as placing con - 
straints on the mass and age of the star dMazumdar et alj |2006). 
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Most recently, iBriquet et al.l d2007l) have successfully mod- 
eled the f3 Cephei star 9 Ophiuchi using rigid rotation. This 
star has bee n the subject of both photometric and sp ectroscopic 
campaigns dHandler et alj|2005t Briquet et al.ll2005l) . with 7 fre- 
quencies identified. These frequencies are thought to be the ra- 
dial fundamental, one ( = 1 triplet, and three components of 
an £ — 2 quintuplet based on both spectroscopic and photo- 
metric mode identification. Recent independent modeling agrees 
with the t identifications, although in some cases they assign 
a diff erent m to the modes dDaszvnska-Daszkiewicz & Walczak 
2010). Spectroscopic observations were also used to determine 
the metallicity of 6 Ophiuchi, w ith a best value of Z = 0.0114 
using the new Asplund mixture (Briq uet et al.ll2007l) . Based on 
these observations, the best fitting model for 9 Oph has been 
found to have a mass of about 8.2 M , X e = 0.38, T e ff = 22 
260, log L/L - 3.85 a nd a rotational velocity of about 30 km 
s" 1 ( IBriquet et al.ll2007l) . 
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The convective core overshoot determined for 9 Ophiuchi is 
surprisingly large, a ov = 0.44 ± 0.07. This result is more than 
double that found for similar /3 Cephei stars, which have a ov 
around 0.1-0.2. It is possible that the unusually high overshoot 
is a result of ID modeling, which does not take into account ro- 
tational effects on the evolution. Although 9 Oph is relatively 
slowly rotating for a B star, its rotation velocity , vsini ~ 30 km 
s _1 , is significantly higher than the stars discussed above, which 
have rotation velocities around 2 km s _1 (V836 Cen) and 6 km 
s _1 (v Eri). Rotation and convective core overshoot may be com- 
plimentary effects, in which case including rotation could reduce 
the amount of convective core overshoot required to match the 
observed frequencies. The vsini of 9 Oph is about 30 km s^and 
the st ar is thought to be viewed nearly equator on (Briquet et al. 
120051) . Although this is not particularly rapid rotation for a star 
of this type, the rotation should produce some effect on the struc- 
ture and frequencies. We have chosen to model 9 Oph using 
models which are uniformly rotating on the ZAMS. In this work, 
we use 2D stellar evolution and linear adiabatic pulsation codes 
to determine the effects of rotation and overshoot on pulsation 
frequencies. By including a 2D treatment of rotation, we investi- 
gate whether the overshoot of 8 Oph could be reduced while still 
matching the observed frequencies. 

This paper is organized as follows. In section [2] we discuss 
the stellar evolution and pulsation calculations. In section[3]we 
recall the effects of convective overshoot on the structure of the 
star, and in section|4]we discuss the effects of rotation and over- 
shoot on the observed frequencies in these models. Using the 
resulting variation we attempt to determine the rotation rate and 
overshoot using the observed frequencies in section [5] As aster- 
oseismic modeling is more commonly done using ID models, 
we perform a "hare and hound" exercise to determine how dif- 
ferent the results from the two methods can be. This exercise 
is discussed in section [6] Finally, we constrain the rotation and 
overshoot of 9 Oph in section Q Our results are summarized in 
section [8] 

2. Numerical Method 

Our stellar m odels are calculated u sing the 2D stellar structure 
code R0T0RC (lDeupreell 990. 1995b. This code takes the conser- 
vation equations for mass, energy, 3 components of momentum 
and the composition together with Poisson's equation and solves 
them implicitly on a two dimensional fin ite difference grid using 



the Henyey method (Henyey et al.lll964l) . We use OPAL opaci- 
ties ([i glesias & Rogers! 19961) andequation of state (Rogers et alJ 



1996). Unlike standard ID stellar evolution codes, ROTORC uses 
the fractional surface equatorial radius and the colatitude as in- 
dependent variables. The surface equatorial radius is determined 
by requiring that the integral of the density over the volume 
of the model equal the total mass. The surface radius in the 
other angular zones is calculated by assuming the surface is an 
equipotential. As the models evolve, angular momentum is con- 
served locally. Although these models are uniformly rotating on 
the ZAMS, they become slightly differentially rotating as they 
evolve. This 2D modeling allows us to calculate the structural 
changes produced by rigid rotation without making assumptions 
about the shape a priori. 

Based on the calculations of Bri quet et al.1 d2007l) . the mass 
of 9 Oph is around 8.5 M . We have found that when rotation is 
taken into account, more massive models are generally needed 
to reach a given temperature and lumniosity. We have calculated 
models at 8.5, 9 and 9.5 M G , and have found that the 9.5 M 
models gave the best match to the observed temperature and lu- 
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Fig. 1. Location in the HR diagram for the 9.5 M models at 
age =15.6 Myr, y eq = (star), 50 (circle), 100 (X), 150 (square) 



and 200 (dot) km s" 



: 0-0.38, increasing from right to left. 



A sample evolution track is shown for a model with v eq - and 
a ov = 0.18. The location and error bars fo r 9 Oph are also show n, 
based on the photometric observations of Handl eret al.l J2005). 



minosity of 9 Oph. We have calculated a small grid of these 9.5 
M models evolved to an age of 15.6 Myr. For these calculations, 
we used 581 radial zones and 10 angular zones. The fractional 
surface equatorial radius was taken to be 1 , and the ratio between 
the polar and equatorial radii decreases as the rotation rate in- 
creases. For the velocities considere d here, the rotational effects 
are not expected to be too severe (see lLovekin et al.l2009t) . These 
models are uniformly rotating on the ZAMS, with surface equa- 
torial velocities on the ZAMS of 0, 50, 100, 150 and 200 km 
s _1 , (corresponding to approximately 0, 35, 65, 1 10 and 150 km 
s at 15.6 Myr) and overshoot parameters (a ov ) of 0, 0.08, 0.18, 
0.28 and 0.3 8. We used a metallic ity slightly higher than that de- 
termined bv lBriquet etail (120071) for 9 Oph, with Z = 0.02. The 
initial hydrogen fraction of these models was 0.7. The proper- 
ties of these models are summarized in Table [TJ The location of 
these models in the HR diagram is shown in Figure[TJ along with 
a sample evolution track (v = 0, a m = 0.18). The location of 9 
Oph with observational error bars from the photometric obser- 
vations of Han dler et al.l (120051) is shown for reference. 

To calculate pulsation frequencies, we used a 2D linear adi- 
abatic pulsation code, NR0 (IClemen afl998h . This code solves the 
linearized pulsation equations on a 2D grid using a finite differ- 
ence technique. The ROTORC model, which is defined on a spher- 
ical polar grid is transformed into a model with the same num- 
ber of radial zones, but defined on surfaces of constant density. 
The pulsation equations are rewritten as finite difference expres- 
sions and the coefficients are placed in a band diagonal matrix. 
Each element of this matrix is itself a matrix, containing the co- 
efficients at each zone in the 2D grid. The solution proceeds in 
two steps, from the centre outwards and from the surface in- 
wards. The solutions are required to match at some intermediate 
fitting surface. At this point, a discriminant can be evaluated, 
which will only be satisfied (equal to zero) if an eigenvalue has 
been located. Frequencies are detected by stepping through fre- 
quency space and looking for zero crossings in the value of the 
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discriminant. This method can result in missed frequencies if 
the frequencies are sufficiently close together, but this can usu- 
ally be avoided by reducing the frequency step size. For furthe r 
discussion of the solution mechanism, refer to IClemenll #1998): 
lLovekin & Deupred d2008l) . 

In rotating stars, a given mode cannot be described by a 
single spherical harmonic, but requires a linear combination of 
Y'"'s. We can calculate the first few terms in this combination as 
the eigenfunctions calculated with NRO are given as a function 
of r and 0, and are defined at several points on the surface. Up to 
9 angular zones can be included in the calculation, with one ra- 
dial integration performed for each angle included. The solution 
is given at N angles, which can subsequently be decomposed 
into the contributions of individual spherical harmonics, effec- 
tively calculating the first N terms in the linear combination. 
Each radial integration contains angular derivatives, evaluated 
using finite differencing, so the resultant coupling among spher- 
ical harmonics arises naturally. In NRO, specifying I specifies the 
parity of the mode, and the calculation is performed using the 
first N even or odd spherical harmonics. In non-axisymmetric 
modes, the appropriate spherical harmonics are chosen starting 
with { — m. 

Since a given mode in a rotating star is described by a linear 
combination of spherical harmonics, the mode no longer has a 
unique I, although m does remain unique. Since t is no longer 
unique, some new way of identifying modes must be used. We 
have chosen to identify modes by the parameter £ a , which is the 
I of the mode in the non-rotating model to which a given mode 
can be traced back. This tracing is done for a series of models of 
gradually increasing velocity, based on both the shape of the lat- 
itudinal variation at the surface and the frequency. As the veloc- 
ity gradually increases, the resulting distortion increases, allow- 
ing the modes to be identified. However, this process becomes 
more difficult as the rotation rate increases, and is described 



in more detail in lLovekin & Deupreel d2008l). Similar pro blems 
have been encountered by for examplejkeese et al. (2009). 

We have calculated frequencies for ( < 2, as shown in 
Figure [2] for the non-rotating models. We have chosen these 
low order p-modes for comparison with the modes detected in 6 
Oph. The frequencies shown in this plot were scaled by a fac- 
tor of -\fiGM/R(40°) 3 ) to account for differences in the radii 
of the models. We have chosen the radius at a colatitude of 
40° as this radius has been found to be the radius most ap- 
propriate for calculating a pulsation constant in similar models 
(Lovekin & Deupredl20 08). 



3. Structural Effects 

In order to better understand the influence of convective core 
overshooting, we have considered the effects of rotation and 
overshoot on the stellar structure for a collection of models at 
the same age. As seen in Figure Q] and Table Q] increasing the 
convective core overshoot increases both the effective temper- 
ature and luminosity slightly. Increased rotation does the op- 
posite, causing a slight decrease in temperature and luminos- 
ity. This decrease is small except for the most rapidly rotating 
model considered here (200 km s~'). These changes in tem- 
perature and luminosity can be related to the change in the ra- 
dius of the star. As the rotation rate increases, the equatorial ra- 
dius increases, while the polar radius decreases. This effect is 
well known, and has been demonstrated by many other authors 
(see, for exam ple iBodenheimerll 1 97 ll iMevnet & Maederlll997l: 
lDeupred20 01). We have also found that for a given rotation rate, 
increasing the convective core overshoot decreases the radius by 
up to 0.8 R Q for models at the same age and velocity, despite the 
corresponding increase in both temperature and luminosity (see 
Table Q]) 
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Fig. 2. The calculated frequencies for the non-rotating model; * 
- € a = 0, o - t Q = 1, x - { = 2. All frequencies are normalized by 
V(GA^/^(40°) 3 ), to account for the differences in stellar radius. 



To help us assess the effect of changing core overshoot on 
the structure, and hence the frequencies, we have calculated the 
Brunt-Vaisala frequency, defined as 



N 2 =. 



1 dP 

pf\!fr 



1 dp 

p dr 



(1) 



for each of our models. In Figure [3] we show N 2 /g for a non- 
rotating model and a rapidly rotating model (200 km s~') with 
no convective core overshooting. The small wiggles visible, par- 
ticularly in the outer sections of the star, are purely numerical 
effects, resulting from the finite difference calculation of the 
derivative. The main peak at the boundary of the convective 
core has approximately the same size and shape in both mod- 
els. Clearly, changing the rotation rate does not produce much 
change in the shape of the normalized Brunt-Vaisala frequency. 
Although increased rotation increases the absolute size of the 
convective core, in terms of fractional radius, this size remains 
roughly constant. In these models, local conservation of momen- 
tum means there is little radial mixing, so the composition gra- 
dient does not change significantly with rotation rate. Although 
there are slight differences throughout the envelope of the star, 
they are small and unlikely to cause large shifts in the frequen- 
cies. On the other hand, the overshoot has a much larger effect on 
the Brunt-Vaisala frequency. Figure [4] shows the Brunt-Vaisala 
frequency for two non-rotating models with core overshoot pa- 
rameters of and 0.38. As expected, the peak at the boundary of 
the convective core has shifted outwards in radius. This region 
is really the only significant difference, and the Brunt-Vaisala 
frequencies are similar throughout the envelopes of the models. 



10" 



-11 



C\2 

S5 



10" 



-12 




0.2 0.4 0.6 0.8 

Fractional Radius 

Fig. 3. The normalized Brunt-Vaisala frequency for a 9.5 M Q 
non-rotating model (solid) and a model rotating uniformly at 200 
km s _1 (dashed). Both models have a convective core overshoot 
parameter of and have been evolved to an age of 15.6 Myr. 
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Fig. 4. The normalized Brunt-Vaisala frequency for 9.5 M non- 
rotating models with core overshoot parameters of (solid) and 
0.38 (dashed), evolved to an age of 15.6Myr. 



4. Effects on the Frequencies and Frequency 
Separations 

Using the 2D stellar evolution and pulsation calculations de- 
scribed in Section [2] we have calculated pulsation frequencies 
for models with rotational velocities from 0-200 km s _1 and 
overshoot parameters of 0-0.38. Each model was evolved to an 



age of 15.6 Myr. As the overshoot and rotation change the evo- 
lution slightly, all of these models have slightly different radii 
and core hydrogen fraction (X c ), as discussed in section [3] (see 
also Table Q] and Figure Q]). In order to offset the effects of 
this difference, we have scaled the frequencies by a factor of 
~\/(GM/R(40°) 3 ) when looking for the effects of rotation and con- 
vective core overshoot. We have chosen the radius at this colati- 
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Fig. 5. The change in the {„ — normalized frequency relative 
to the frequency of the non-rotating model (stars) as a function 
of overshoot for models rotating at 50 km s _1 (circles), 100 km 
s _1 (X) 150 km s _1 (squares) and 200 km s _1 (triangles) for the pi 
(bottom) and p2 (top) modes. 
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Fig. 6. The change in the l — normalized frequency relative 
to the a ov = case (stars) as a function of velocity for the pi 
(bottom) and p2 (top) modes. Shown are overshoots of a ov = 
0.08 (circles), a ov - 0.18 (X), a ov = 0.28 (squares) and a ov - 
0.38 (triangles). 



tude as it has been found t o produce the best pulsation constant 
(lLovekin & De upree 2008), and hence should be the best choice 
for scaling the frequencies. 

4.1. Frequencies 

To determine the change in frequency produced by rotation and 
overshoot, we have plotted the scaled frequency as a function of 
overshoot relative to the non-rotating case (Figure [5]) for the t 
— pi and p2 modes. The results are similar for the t = \ and 2 
modes. In all cases, the differences in scaled frequency produced 
by increasing the velocity are small, of order 1-1.5%. This is true 
regardless of which radius is used to scale the frequencies. Using 
the equatorial radius does not change the results, while when the 
frequencies are scaled by the polar radius the magnitude of the 
frequency difference is slightly larger. Interestingly, this change 
in scaling also changes the sign of the frequency differences. 
When the frequencies are scaled by the equatorial radius, the 
frequencies increase as the rotation rate increases, while they 
decrease with increasing rotation rate when scaled by the polar 
radius. 

As discussed in section [3] although the equatorial radius in- 
creases as the rotation rate increases, the polar radius can actu- 
ally decrease slightly. For example, while the equatorial radius 
in the models with no convective overshoot increases from 6.76 
R to 6.85 R Q as the rotation rate increases from to 200 km s _1 , 
the polar radius actually decreases to 6.68 R . This increases the 
scaling factor as the rotation rate increases, and so the scaled fre- 
quencies actually decrease relative to the non-rotating case. At 
a given overshoot, the radius at a colatitude of 40° stays nearly 
constant, changing by only 0.007 R Q as the velocity increases 
from to 200 km s" 1 , giving us another reason to choose this as 
the scaling radius. However, for a given velocity there is no point 
on the stellar surface at which the radius stays approximately 
constant for all overshoots, and scaling by the radius is not as 
effective. In this way, our scaling accounts for the change in ra- 



dius produced by rotation, but not that produced by convective 
core overshooting. As a result, the curves shown in Figure|6]are 
relatively flat as a function of velocity, while the curves shown in 
Figure show variation with increasing overshoot. This choice 
of scaling will highlight the differences resulting from the over- 
shooting, making them easier to analyze. 

We have also calculated the change in frequency as a func- 
tion of velocity relative to the case with a ov = 0, shown in Figure 
[6] again for the l — pi and p2 modes. The trends are similar, 
although larger in magnitude, with the differences typically of 
order 10 %, as might be expected based on our choice of scaling. 
In this case, overshoot appears to affect all t values equally, as 
there is little variation with increasing t . When scaling the fre- 
quencies, using the polar radius decreases the magnitude of the 
frequency difference slightly, but does not affect the sign of the 
frequency difference. Changing the convective core overshoot 
does not affect the shape of the star as dramatically as does ro- 
tation, and both the polar and equatorial radii decrease as a ov is 
increased. 

Both velocity and overshoot can significantly affect the ab- 
solute frequencies, although the effect is almost negligible (less 
than 1.5%) with respect to rotation when the frequencies are 
scaled to account for differences in radius. The differences 
shown in Figures [5] and [6] are scaled using the radius at 40°; 
the differences in the unsealed frequencies are about two orders 
of magnitude larger. Although the scaled frequencies decrease 
with increasing rotation and overshoot, increasing overshoot ac- 
tually increases the unsealed frequencies, while increasing ve- 
locity causes the frequencies to decrease. Based on the results 
here, we find that the effects of slow to moderate uniform ro- 
tation (up to about 200 km s _I ) can be primarily accounted for 
by the effect on the stellar radius. The effects of increasing con- 
vective core overshoot, on the other hand, are not so easily ac- 
counted for by scaling, and are nearly an order of magnitude 
larger than the rotational differences when scaled frequencies 
are considered. Although the value of the frequencies themselves 
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Fig. 7. The large separation for the l — 2 pi and p2 modes as 
a function of overshoot for velocities of km s _1 (stars), 50 km 



Fig. 8. The small separations for even (l = 0, 2, p3 and P4) 
modes as a function of overshoot for km s _1 (stars), 50 km 
s _1 (circles), 100 km s _1 (X), 150 km s _1 (squares) and 200 km s _1 (circles), 100 km s _1 (X), 150 km s _1 (squares) and 200 km 
s" 1 (triangles). All frequencies are scaled by y[{GMIR y ) to ac- s (triangles). As for the large separation, the frequencies are 



count for differences in radius. 

can change significantly, we find that as might be expected, over- 
shoot produces no effect on the mode splitting. 



4.2. Frequency Separations 

We have also calculated the effect of rotation and overshoot on 
the frequency separations. Although the large and small sepa- 
rations are asymptotic limits, we have applied the definitions to 
our low order modes, and continue to refer to these as the large 
and small separation. At large radial order, the large separation 
is defined as 



Av, 



l - V{,n+l - Vtji, 



(2) 



where n is the radial order of the modes. This separation is shown 
in Figure|7]for the low order ( =2 modes. Both increasing ro- 
tation and increasing core overshoot produce a decrease in the 
large separation, most dramatically for the t = 2 modes. As the 
overshoot increases, the large separation decreases by about 7- 
10 % for the l = and 1 modes, but falls by nearly 50 % for the 
I — 2 modes. The decrea se in large separation with increasing 
rotation rate, as found bv lLovekin et alJ (120091) is also seen, al- 
though the effects are not as dramatic as the effects of convective 
core overshoot. Again, this is at least in part due to our choice of 
scaling, which minimizes the effects of rotation to better show 
the effects of core overshoot. 

The small separation is defined as 



d(,n - V(,n - Vf+2,n-l 



(3) 



and shown in Figure [8] for the even modes. As with the large 
separation, the small separation decreases slightly as a function 
of rotation rate up to approximately 200 km s . At this veloc- 
ity, there is a minimum in the small separation, and the small 
separation increases by appro ximately 20 yuHz as the rotation in- 
creases beyond this velocity (lLovekin et alj2 009). The effects of 
overshoot on the small separation are shown in Figure [8] were a 
significant increase can be seen as the core overshoot increases. 



scaled by ^fiGM/P?) to account for the different radii of the 
models. 



A large increase in the small separation could therefore be in- 
dicative of either a large core overshoot or a high rotation rate. 



5. Constraining a ov 

Once we had established the effects of rotation and overshooting 
on the pulsation frequencies and separations, we attempted to 
determine if these could be used to constrain a star given a set of 
observed frequencies. To do this, we calculated three test models 
with the same mass and age as the previous grid of models (see 
Section but with slightly different rotation rates and convec- 
tive core overshoots: (y eq ,a in ) = (80 km s _1 , 0.1), (180 km s _1 , 
0.1) and (80 km s _1 , 0.3). By comparing the differences between 
the frequencies for these test models (henceforth the observed 
frequencies) and the frequencies in our original grid (model fre- 
quencies) we hope to constrain the rotation and overshoot. Note 
that throughout this section, we compare absolute frequencies, 
not scaled frequencies as in the previous section. 

Initially, we attempted to match individual frequencies. For a 
given observed frequency, we calculate the absolute value of the 
difference between the observed frequency and each model fre- 
quency {\vt, g ridlyt=o^bs ~ Ve,ob s /V{=o,ob s \) and plot this difference 
versus the observed frequency (ye <o bslvt=0,obs)- All frequencies 
are normalized by the observed i„ — frequency to offset the ef- 
fects of varying radius. The results of this calculation are shown 
in Figure|9]for the l — pi frequency of the observed M eq = 80 



km s 



0. 1 case. The models with the same convective core 



overshoot parameter cluster in this diagram, giving no informa- 
tion about the velocity. There is a clear minimum difference in 
frequency which falls between a ov - 0.08 and 0.18, correspond- 
ing to the target overshoot of 0. 1 . 

Using this method, we can of course find a good match for 
a single frequency. However, we found that fitting different ob- 
served frequencies from a single observed star can give different 
results for the best fitting overshoot. To determine the single best 
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Fig. 9. Determination of overshoot using the ( — pi frequency 
for an observed model with a = 0. 1 , v = 80 km s~ 1 . The observed 
frequency is compared to the corresponding frequency for each 
model in the grid, plotted against the observed frequency. All 
frequencies are normalized by the observed („ = frequency. 
The models corresponding to constant a ov cluster together in this 
diagram: a ov = (stars), a ov = 0.08 (circles), a ov = 0.18 (X), a ov 
= 0.28 (squares) and a ov =0.38 (stars). The best fit model has 
a ov between 0.08 and 0.18. 



solution for the rotation and overshoot for each set of observa- 
tions, we combine the data for all observed frequencies to deter- 
mine the best overall fit for a given model. We do this using a 
modified^ 2 statistic. To calculate this, we compare the observed 
frequencies to the calculated frequencies from each model. For a 
given rotation and overshoot, we determine the best fitting model 
frequency for each observed frequency, where the ( a of the mode 
is known, but m and n are free parameters. The total x 2 for this 
model is the sum of the differences squared for the best-fitting 
matches for the individual frequencies, normalized by the num- 
ber of frequencies we fit: 



X 



(4) 



where N is the number of frequencies fit for each model. The sur- 
face equatorial velocity and core overshoot of the best fit model 
is given by the lowest overall x 2 - The x 2 calculated for each 
model can then be plotted as a function of velocity and over- 
shoot, as shown in Figure [10] giving us a simple visual way of 
determining the best fitting model. 

This method works best if some additional constraints are 
applied. Once a model frequency has been matched to an ob- 
served frequency, it is flagged and considered unavailable for 
future matches. If more than one observed frequency is best 
matched to the same model frequency, we use a recursive algo- 
rithm to determine which arrangement of the frequencies gives 
the lowest x 2 ■ This also allows us to determine the best fit in 
cases where we have more observed frequencies than model fre- 
quencies. We also exclude observed frequencies that are more 
than 5 /jHz above (below) the maximum (minimum) value in the 
model grid. Both of these constraints remove frequencies from 
the set of observed frequencies, and different numbers of fre- 
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Fig. 10. The log of^ 2 for an "observed" model with a ov = 0.3 
and v eq = 59 km s _1 . The minimum can be seen around (65,0.28), 
with models at higher velocity and similar overshoot also ac- 
ceptable. Although the resolution is course, we can place some 
constraints on a given model. 



quencies may be removed when comparing to different models. 
To account for this, for each overshoot and velocity, we nor- 
malize the x 2 by the number of observed frequencies for which 
matches were found. 

Using this method, we were able to recover the input pa- 
rameters for our three "observed" models. The observed models 
have velocities of 80 and 180 km s _1 on the ZAMS, which corre- 
spond to 57 or 65 km s _1 (depending on the core overshooting) 
and 129 km s~'at an age of 15.6 Myr, and overshoot parame- 
ters of 0.1 and 0.3 H p . For a model with (v = 129 km s _1 , a ov 
= O.lHp), we recovered (110, 0.08), with an uncertainty of + 50 
km s _1 in velocity and ± 0.1 in overshoot due to the coarseness 
of our model grid. Similarly, for a model with (57, 0.1) we re- 
covered (65, 0.18) and for a model with (59, 0.3), we recovered 
(65, 0.28). A sample plot of the^ 2 space is shown in Figure [TOl 
for a model with parameters (59, 0.3). In general, we find accept- 
able solutions over a range of bins, generally straddling the true 
value. In the example shown in Figure[l0] the best solution has a 
slightly higher velocity and lower overshoot than the true values, 
and models with a slightly lower velocity or higher overshoot are 
also acceptable solutions. 

The models shown above were calculated using a full set of 
calculated frequencies (about 20 frequencies), including the first 
two radial orders for all modes. In more realistic cases, fewer 
frequencies would be available, and the m and n of these modes 
would be unknown. To examine a more realistic case, we have 
created an "observed" star with a reduced set of frequencies and 
compared them to the full model grid. We have assumed the fre- 
quencies of the I — fundamental, a t — 1 triplet and a full 
I — 2 quintuplet are known (10 frequencies in total). We found 
that reducing the number of frequencies in this way did not af- 
fect the determination of the best fitting solution. We do find 
some evidence for a degeneracy between rotation and overshoot 
in the most rapidly rotating case, shown in Figure [TT] In this 
case, when only 10 frequencies are included, the best fitting so- 
lution remains in the same area, now at 145 km s _1 , a ov = 0.08. 
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Fig. 11. The log of x 1 for an "observed" model with a ov = 0.1 
and v eq = 129 km s _1 , fit using 10 frequencies. The best fitting 
solution can be seen at (145 km s _1 , 



0.08), while a secondary 
minimum has appeared at lower velocity and higher overshoot. 

A secondary minimum has also appeared, with acceptable solu- 
tions having a core overshoot parameter of 0.38 and a velocity 
of 65 km s _1 . 

Clearly, our resolution, particularly in velocity, is very low, 
and better constraints could be established by using a finer grid. 
However, even with this coarse resolution, the method success- 
fully returns the correct (albeit uncertain) surface equatorial ve- 
locity and core overshoot, at least in cases where the mass and 
age of the star are known. In general, the mass and age are not 
known, which will make the fitting process more difficult. 

6. Hare and Hound Exercise 

Next, we decided to see how well this determination of rotation 
and overshoot could be done using ID models. This will help 
quantify the errors introduced by neglecting rapid rotation in 
these models. To do this, low order frequencies were calculated 
by one of us (CCL) using a uniformly rotating 2D model as de- 
scribed above. The selected model is an 8.5 M model, evolved 
to an age of 20 Myr, with a ZAMS rotation rate of 150 km s _1 
(108 km s _1 after 20 Myr). The calculated frequencies and ob- 
served position in the HR diagram were passed on as a set of 
observational quantities, given in Table|2] One of us (MJG) then 
used these observations to try and determine the convective core 
overshoot and rotation rate of the star. The star was assumed to 
be a uniformly rotating [3 Cephei star oscillating with low order p 
and g modes with metallicity Z = 0.02. For comparison, we also 
use the frequencies in Table|2]in a 2D fitting process as described 
above, assuming the mass and age of the star are known. 

In the ID case, this star was modelled with masses of 8.3 and 
8.5 M , chosen based on which tracks cross the observed loca- 
tion in the HR diagram. The calculated evolutionary sequences 
cross the error box at an age between 18 and 25 Myr, depend- 
ing on the amount of convective core overshoot. Sequences were 
considered for a — 0, 0.05, 0.1, 0.2, 0.3 and 0.4. As expected, 
there is a degeneracy between mass, age and overshoot in the HR 
diagram, shown in Figure [12] The oscillation code used in this 
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Fig. 12. HR diagram showing the evolutionary tracks (dashed 
lines) for 8.5 M (black) and 8.3 M (red) models with different 
values of the core overshoot parameter a ov . Open dots represent 
models at time intervals of 1 Myr around an age of 20 Myr (filled 
dots). The box represents the uncertainty of the observed loca- 
tion of the "star" in this diagram, as given in Table [2] 



section, WAR(saw)M(eudon), includes the effects of uniform ro- 
tation as a perturbation for axisymmetric modes as well as fre- 
quencies of nonaxisymmetric modes up to third order in the rota- 
tion rate ( Soufi et alJll998UDaszvnska-Daszkiewicz et ai1 [2002; 
iGoupil & Talod2009tlGoupii2009l.l20ldr Comparison between 
frequencies obtained with a pertur bed approach and a nonpertur- 
bative one can be found in Reese et all J2006): Lignieres et al.l 



(2006) for a frequency comparison usi ng the same model, a 2D 
poly trope, in both approaches, and in lOuazzani et alJ d2009l) for 
a comparison using a ID poly trope in the perturbative approach 
and a 2D poly trope involved in the nonperturbative one. 

Initially, frequencies were calculated assuming no rotation. 
Particularly for the i = 1 modes, matching the frequencies was 
a problem, and for some modes no solution was found. In this 
case, the modes identified in Table [2] as pi, pi and g\ had fre- 
quencies corresponding to the p\, gi and g2 modes in the ID 
models. None of the models considered were able to simulta- 
neously fit all of the modes at the same age. A x 2 technique, 
similar to that described above, was then used to fit all of the 
modes simultaneously. This x 2 fit determined that the 8.3 M 
models gave a better fit to the data, with overshoots of around 
0.1 ± 0.05. The models with the lowest^ 2 have ages between 
21-22 Myr, and the temperature and luminosity are within the 
observed errors of the target "star" (see Figure [T3l. 

The ID models used in this section are capable of including 
the effects of uniform rotation as a perturbation, allowing the 
mode splitting to be included. A second x 2 was calculated for 
the splitting of the models compared to the "observed" model. 
As expected, this splitting was found to be independent of the 
overshoot. Good agreement was found for rotation rates below 
125 km s _1 , while above this rate the g mode splitting does not 
agree and^- 2 becomes large. 
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Fig. 13. The quantity^ 2 for modes {( — p\, pi, and { = 1 g\, 
g2, and pi) as a function of overshoot for models with mass of 
8.5 M (blue) and 8.3 M (red). 



When all effects are included, the ID modeling finds the best 
solution to be 8.3 M Q , 21 Myr, a ov =0.1 H p and \ eq = 125 km 
s _1 . In contrast, the 2D modelling described above finds a best fit 
model with a ov = 0.28 and v eq = 145 km s _1 , given an assumed 
mass and age of 8.5 M Q and 20 Myr. Solutions at 65 km s _1 are 
also acceptable, as are core overshoot parameters of 0.18. The 
large range in possible velocities found with the 2D models are 
a result of the coarseness of the grid, which contains no mod- 
els with velocities near 100 km s _1 . Although the ID models 
are able to determine the velocity of a target model, the deter- 
mination of the oversh oot i s much less accurate . As di scussed in 
llignieres etail (120061) and lLovekin & Deupred (120081) . for stars 
rotating more rapidly than about 100 km s , a more accurate 
treatment of the rotation and pulsation is needed. These works 
have shown that a 2D approach tends to produce a larger effect 
on the frequencies than ID perturbative methods. As discussed 
above, rotation and overshoot change the unsealed frequencies 
in opposite directions. While increasing rotation decreases the 
frequencies, increa sing core overshoot increases the frequen- 
cies. A ccording to Lig nieres et al.l d2006l) : lLovekin & Deu pree 
(120081) . ID calculations will not find as large a shift from rota- 



tion as 2D models, and hence will need a smaller overshoot to 
match the observed frequency in rapidly rotating stars. 

7. Application to 9 Oph 

Finally, we decided to test the technique discussed in Section 
[5] on 9 Oph. As discussed in Section [5] our resolution in both 
velocity and overshoot is quite coarse, so we do not expect to 
place tight constraints on 9 Oph in this exploratory calculation. 
We fixed the metal l icity o f our models at Z = 0.02, using the 
iGrevesse & Sauval d 19981) abundances. For our asteroseismic 
comparison, we fixed the mass at 9.5 M Q , after calculating the 
evolutionary tracks of 8.5, 9 and 9.5 M models and determining 
which models gave the best match to the observed location of 9 
Oph. 

The seven observed frequencies of 9 Oph were compared 
against the 9.5 M models shown in Figure [1] evolved to an 
age of 15 .6 Myr. We used the spectroscopic mode identifications 
given by Briq uet et al.l ([2005), but used only the given £ values 
as constraints. As for the calculations performed in Section|5] we 
found that the best results were obtained when each frequency is 
compared only to model frequencies of the same {„. When we 
do this comparison for 9 Oph, we find the best fit model has a 
rotation and overshoot of (35 ± 35 km s _1 , 0.28 ± 0.1H p ). This 
model is located within the observed photometric errors for 9 
Oph, with a luminosity of log(L/L Q ) = 3.75 and an effective 
temperature T e ff = 225907^. The rotation velocity found in our 
best solution is in good agreement with the rotation velocity as 
determined fro m both the vsini and the o bserved mode splitting, 
29 ± 7 km s _1 dBriquet et al.ll2005l 120071) . The convective core 
ov ershoot determined h ere is slightly lower than that determined 
by Brique fet al.l (120071) . although the two results do agree within 
the errors. It seems that when a more accurate treatment of rota- 
tion is taken into account, at least for slowly rotating stars, this 
does indeed reduce the need for convective core overshooting. 
Our new estimate is still higher than the results determined for 
other [3 Cephei stars (typically 0.1-0.2), although again, the re- 
sults agree to within our uncertainties. 

Our models for 9 Oph have been calcu lated using a metallic - 
ity of Z = 0.02 at the composition of lGrevesse & Sauval (fl998h 
consid erably higher than that determined by Briquet et al. 
(120071) . As discussed in their work, higher metallicity is corre- 
lated to lower convective core overshoot, so one would expect 
our models to have a lower overshoot. In fact, based on the re- 
sults given in their Table 4, we would expect a ID model with our 
metallicity to have an overshoot of a ov = 0.32. A s is expected, 
this is lower than the overshoot calculated by Brique t~et al.l 
(2007]), and is much closer to the value obtained with our 2D 
models. Nevertheless, our models also show the need for an un- 
usually high convective core overshoot in this star. Although our 
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Overshoot 

Fig. 14. Log x 2 plotted as a function of velocity (y-axis) and 
convective overshoot (x-axis). The minimum (green) in x 1 for 
this star is for a = 0.28, with a velocity of 35 km s . 

models have shown some indications that rotation and overshoot 
can be complimentary effects, it appears that more rapid rotation 
is required in order to influence the frequencies. 

Another limitation of our models is the restriction to fixed 
mass and age. This is of course not realistic, and has been 
introduced to limit the computations required. As mentioned 
in Section [2] we initially calculated models at three different 
masses, 8.5, 9 and 9.5 M . All of these models were close to 
the observed error box of 9 Oph, although the lower mass mod- 
els needed to be evolved to a higher age. Although detailed 
frequency calculations were not performed for these models, it 
seems likely that matching models could be found. Including 
mass and age as free parameters is currently under investigation 
and will be included in future work. 

8. Conclusions 

We have found that both rotation and overshoot do have an effect 
on stellar frequencies, although this is predominately through 
the effects of changing stellar radius. Increasing convective core 
overshoot increases the unsealed frequencies, while rotation 
causes the unsealed frequencies to decrease, although certain 
choices of scaling can minimize this effect for rotation. Although 
the frequencies themselves may change, as expected we find that 
the overshoot has no effect on the mode splitting. 

We also investigated the effect of overshooting on the large 
and small separations. We find that increasing convective core 
overshoot decreases the large separations, but increases the small 
separation. Unfortunately, bo th of these effects ar e similar to the 
effects of rotation noted by lLovekin et al.l d2009h . and may not 
be useful for disentangling the effects of rotation and overshoot. 

We have attempted to use the changes in frequency to con- 
strain the best fitting convective core overshoot. Using individual 
modes, we were able to easily constrain the convective core over- 
shoot, although not the rotation. As should be expected, differ- 
ent modes gave different results, so we developed a modified ;f 2 
statistic, simultaneously fitting all known frequencies for models 
of known mass and age (9.5 M G , 15.6 Myr). We found this was 



most effective if the i s of all frequencies were known. When 
the £ s are included, we are able to correctly determine both ro- 
tation and overshoot for models in our grid as well as for three 
test models with slightly different v eq and a ov . Although the un- 
certainties on our results are large as a result of the coarse grid 
spacing used, in all cases the best solutions surrounded the true 
parameters of the test models. 

We also conducted a hare and hound exercise using frequen- 
cies calculated using 2D stellar models and pulsation calcula- 
tions. We found that using ID stellar models, we were able 
to find a model that reproduced the frequencies for reasonably 
close values of age, mass, rotation and convective core over- 
shoot. The ID models were able to find a velocity quite close to 
the true rotation rate, although the determination of core over- 
shoot was not as accurate. For a model with true parameters 
(108, 0.25), the ID models returned a best fit of (125, 0. 1), while 
the 2D models found a best fit of (145, 0.28). The mass and age 
returned by the ID models were also quite close to the true val- 
ues. 

Finally, we applied these methods to 9 Oph, a rotating fi 
Cephei star with seven observed frequencies. Using models at 
a fixed mass and metalli city, we find an overs hoot slightly lower 
than that determined by Bri quet et alj d2007l) . with a - 0.28 ± 
0.1, as expected for our higher metallicity models. The rota- 
tion velocity for this model agrees well with the observed value, 
around 30 km s _1 . We also find that a good ma tch to 9 Oph 
requir es a more massive star than determined by Briquet et al. 
(l2007h . around 9.5 M Q vs 8.2 M in their models. This differ- 
ence in mass could be a a consequence of including rotation in 
our models, which tends to make models appear less massive, 
but may also be a result of the higher metallicity used in our 
models. Decreasing the convective core overshoot from 0.44 to 
0.28 H p brings 9 Oph closer to the range of overshoots found in 
other p Cephei stars, around 0. 1-0.2 H p , but is still high. 

As discussed above, (Section[6]), the fact that the ID calcu- 
lation returns a lower overshoot than the 2D models while the 
ID calculations find a higher overshoot for 9 Oph is probably a 
result of the ID treatment of rotation. At slow rotation velocities 
(< 100 km s~'), like 9 Oph, the ID method works well, while 
for more rapidly rotating models, as in the hare and hound ex- 
ercise performed here, the ID treatment of rotation results in an 
underestimate of the convective core overshoot. 
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